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ABSTRACT 

We analyse the stochastic background of gravitational radiation emitted by a cos- 
mological population of core-collapse supernovae. The supernova rate as a function 
£N) ■ of redshift is deduced from an observation-based determination of the star formation 

' rate density evolution. We then restrict our analysis to the range of progenitor masses 

OS ! leading to black hole collapse. In this case, the main features of the gravitational-wave 

■ emission spectra have been shown to be, to some extent, independent of the initial 

conditions and of the equation of state of the collapsing star, and to depend only 
on the black hole mass and angular momentum. We calculate the overall signal pro- 
duced by the ensemble of black-hole collapses throughout the Universe, assuming a flat 
cosmology with vanishing cosmological constant. Within a wide range of parameter 
values, we find that the spectral strain amplitude has a maximum at a few hundred 

■ Hz with an amplitude between 10 -28 and 10 -27 Hz -1 / 2 ; the corresponding closure 

density, SIgWi has a maximum amplitude ranging between 10 -11 and 1CT 10 in the fre- 

q . quency interval ~ 1.5 — 2.5 kHz. Contrary to previous claims, our observation-based 

!— i ' determination leads to a duty cycle of order 0.01, making our stochastic backgound 

, a non-continuous one. Although the amplitude of our background is comparable to 

d • the sensitivity that can be reached by a pair of advanced LIGO detectors, the charac- 

teristic shot-noise structure of the predicted signal might be in principle exploited to 
design specific detection strategies. 

■ Key words: gravitational wave background - star formation rate: black holes. 



1 INTRODUCTION 

According to the theory of general relativity, gravitational 
waves are expected to be produced in such a large variety of 
astrophysical and cosmological phenomena that it is plausi- 
ble to guess that our universe is pervaded by a background 
of gravitational radiation. Depending on their origin, differ- 
ent contributions to this background will exhibit peculiar 
spectral properties, the features of which it is interesting to 
investigate in view of a possible detection by future gravita- 
tional antennas. 

Gravitational waves of cosmological origin could be the 
result of processes that developed in the very early Universe. 
Due to the extremely small graviton cross-section, the spec- 
tral properties of such a relic radiation should have been 
retained with no substantial alteration until today, and, if 
observed, would provide information on the physical condi- 



tions of the epochs when these waves were produced (see, 
for a recent review, Maggiore 1998). 

Further contributions to the stochastic background of 
gravitational waves are of astrophysical nature and their 
generation dates back to more recent epochs, when galax- 
ies and stars started to form and evolve. For instance, the 
background generated by the radiation emitted from rotat- 
ing neutron stars in the Galaxy has been considered in (Gi- 
ampieri 1997; Postnov 1997) and it has been shown that, by 
taking into account the uncertainty with which the elliptic- 
ity of neutron stars is known, /i 2 Ogw(^ = 100 Hz) is in the 
range 10 -15 — 10 -9 , where Q.gw is the closure energy density 
of gravitational waves per logarithmic frequency interval, 
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with p c = 3Hq 1 8nG the critical density.^] A strategy for 
detecting these signals with one interferometric antenna, by 
squaring the detector amplitude and searching for a sidereal 
modulation has been proposed in (Giazotto, Bonazzola & 
Gourgoulhon 1997; Giampieri 1997). Further studies have 
estimated the intensity of the stochastic background pro- 
duced by the galactic merging of unresolved binary white 
dwarfs (Postnov 1997), and subsequently extended to the 
extragalactic white dwarf's merging for different cosmologi- 
cal models (Kosenko & Postnov 1998). The considered range 
of frequency is 10~ 3 — 10~ 2 Hz. 

In this paper, we study the background of gravitational 
waves emitted by sufficiently massive stars that, having 
reached the final stages of their evolution, collapse to form 
a black hole. We consider only collapses to a black hole be- 
cause numerical studies have shown that the spectrum of 
the gravitational energy emitted in these processes exhibits 
some distinctive features that are to some extent indepen- 
dent of the initial conditions and of the equation of state of 
the collapsing star, and that depend only on the black hole 
mass and on its angular momentum (see, for a recent re- 
view, Ferrari & Palomba 1998). Conversely, the simulations 
of gravitational collapse to a neutron star available in the lit- 
erature indicate that the energy spectrum which is produced 
strongly depends on the equation of state of the collapsing 
star, and no distinctive features can easily be extracted to 
model the process. 

The energy spectrum which we will use to model the 
collapse to a black hole is that obtained by Stark & Piran 
(1985, 1986), who studied the evolution of the axisymmetric 
collapse of polytropic, rotating stars, by a fully relativistic 
numerical simulation. 

In order to evaluate the contribution that each burst of 
gravitational radiation emitted in a collapse occurred in re- 
cent or past epochs gives to the stochastic background, the 
knowledge of the star formation rate as a function of the cos- 
mological redshift is essential. This function has been calcu- 
lated in the framework of the hierarchical theory of galaxy 
formation (Cole et al. 1994; Kauffmann, White & Guider- 
doni 1993; Navarro, Frenk & White 1996). However, in this 
paper we prefer to base our calculation on the star forma- 
tion rate which can be deduced from observations, being 
aware of the fact that, since these observations are reliable 
only up to redshifts of order z ~ 4 — 5, we may overlook 
the possibility that a big burst of star formation may have 
occurred in a far past (Rosi & Zimmermann 1976, Bertotti 
& Carr 1980, Bond & Carr 1984). Following a recently pro- 
posed approach (Madau, Pozzetti & Dickinson 1997), the 
star formation rate history that we have considered is based 
on a very simple stellar evolution model which assumes a 
time dependent star formation rate and a constant initial 
mass function. The model is directly inferred from the im- 
pressive collection of recent observations of the rest frame 
UV-optical emission from star forming galaxies at redshifts 
ranging from zero up to z ~ 4 — 5 (Steidel et al. 1996; Madau 
et al. 1996; Lilly et al. 1996; Connolly et al. 1995; Ellis 1997), 
and assumes a Salpeter initial mass function, a flat cosmol- 

* Throughout this paper we shall write the Hubble constant as 
H = hx 100 km s _1 Mpc _1 . 



ogy with zero cosmological constant and h — 0.5. The result- 
ing star formation history has been extensively investigated 
and proves to correctly reproduce, within the experimen- 
tal errors, the evolution of the galaxy emission properties 
observed in other wavebands. Furthermore, it is consistent 
with the indications of QSO absorption lines and metallic 
clouds observations (Ellis 1997; Madau, Pozzetti & Dickin- 
son 1997; Lanzetta, Yahil & Fernandez-Soto 1996; Pei & Fall 
1995). 

The plan of the paper is as follows. In Section 2 we 
briefly review the current state of observations of the global 
star formation in field galaxies and present the model for 
the star formation rate evolution that we have adopted in 
our analysis; in Section 3 we compute the rate of collapses 
and the duty cycle of the process. In Section 4 we describe 
the energy spectrum computed by Stark and Piran, which 
we use to model the gravitational emission of single events. 
In Section 5 we derive the expression for the spectral energy 
density and for Qgw, and discuss the spectral properties 
of this background. In Section 6 we analyze the dependence 
of our results on the assumed model parameters. Finally, in 
Section 7 we present a preliminary statistical analysis of the 
signal. A discussion of the assumptions on which our results 
are based and on further developments of the present work 
will be given in the concluding remarks. 



2 THE STAR FORMATION HISTORY 

In the past few years, our understanding of the origin and 
evolution of galaxies has greatly improved. Within the gen- 
eral framework of hierarchical clustering and gravitational 
instability, numerical and analytical modelling describe how 
primordial perturbations turn into visible stars (Cole et al. 
1994; Kauffmann et al. 1994; Navarro et al. 1996). Notwith- 
standing the complex physical processes at work, the current 
theoretical models succeed in reproducing some of the ob- 
served properties of the galaxy population. However, many 
fundamental aspects still remain poorly understood (Frenk 
et al. 1996; Baugh et al. 1997). 

The major progress made so far comes from the obser- 
vational work. The spectacular new data obtained with the 
Hubble Space Telescope (HST), Keck and other large tele- 
scopes have extended our view of the universe up to z~ 4—5 
(Steidel et al. 1996; Madau et al. 1996; Lilly et al. 1996; 
Connolly et al. 1995; Ellis 1997). 

In the following, we have attempted to summarize the 
relevant aspects of the star formation history that is emerg- 
ing from the data collected so far, having in mind a reader 
interested in gravitational waves but not necessarily compe- 
tent in galaxy evolution. 

The identification of star-forming galaxies at 2 ~ z ~ 4 in 
ground-based surveys and in the Hubble Deep Field (HDF) 
was successfully accomplished through the elaboration of 
robust selection criteria based on multicolor broadband ob- 
servations of the emitter's rest-frame UV and optical stel- 
lar continuum. Ground-based observations of star-forming 
galaxies at z ~ 3 have used color techniques which are sen- 
sitive to the presence of a Lyman-break superimposed on 
an otherwise flat UV spectrum (Steidel & Hamilton 1992). 
At higher redshifts, the effect of intergalactic absorption on 
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galaxy colors can severely affect galaxy spectra. Its inclusion 
in the construction of reliable selection criteria proved ex- 
tremely efficient in identifying star-forming galaxies at red- 
shifts 2 < z < 3.5 (UV-dropout technique) and 3.5 < z < 4.5 
(blue-dropout technique) in the HDF images (Madau et al. 
1996). 

The combination of these high redshift observations 
with the recent completion of several ground-based spec- 
troscopic surveys out to z ~ 1 (Lilly et al. 1996; Ellis et al. 

1996) have enabled, for the first time, a systematic study of 
field galaxies at increasing cosmological lookback times. Fo- 
cusing on the entire population, a great deal can be learned 
about the evolution of galaxies from the analysis of the in- 
tegrated light emitted at each redshift in a given waveband, 
i.e., the comoving luminosity density, 

Px{z) = I L x <t>{L x ,z) dL x , (2) 
Jo 

where <f> (L\, z) is the best-fit Schechter luminosity function 
in each redshift bin. This quantity is independent of many 
details of galaxy evolution, such as the merging history or 
short-lived star formation episodes of individual galaxies, 
and depends on the global star formation history and on 
the initial mass function (IMF) of the stars. In particular, 
the UV continuum emission from an actively star-forming 
galaxy is mainly contributed by massive (M > 10 M©) short- 
lived (tMS < 2 x 10 7 yr) stars. It has been shown that, after an 
initial transient phase, a steady state is reached where the 
measured luminosity becomes proportional to the star for- 
mation rate (SFR) and independent of the galaxy history for 
t 2> tMS (Madau, Pozzetti & Dickinson 1997). Assuming a 
universal IMF, the UV-SFR relation is determined through 
specific stellar population synthesis codes (Bruzual & Char- 
lot 1993; 1998) and the global star formation rate evolution 
is directly inferred from the observed UV emissivity (Madau 
1997; Lilly et al. 1996). 

The comoving SFR density (the mass of gas that goes 
into stars per unit time and comoving volume element) as a 
function of redshift, p*(z), is shown in Figure 1 for a Salpeter 
IMF and a flat cosmology (fi = 1, A = 0), with h = 0.5. 

Despite the uncertainties due to the connection of dis- 
parate datasets, each using a different indicator of the star 
formation activity and each selecting with a specific criterion 
only some detectable subset of the overall population, the 
observations collected so far appear to provide a consistent 
picture of the global star formation history of the universe. 
The SFR density rises sharply from its local value to a peak 
at z ~ 1.5 to fall again out to z as 4. This evolutionary be- 
haviour indicates that the bulk of the stellar population was 
assembled in the redshift range 1 ~ z ~ 2 (Ellis 1997; Madau 

1997) . 

One major source of uncertainty is the effect of dust 
extinction. At redshift z ~ 3, galaxies are already enriched 
in heavy elements and it is likely that some dust is mixed 
with gas and stars. Even a relatively small amount of dust 
can attenuate the UV luminosity and reradiate the absorbed 
UV light in the far-IR thus leading to an underestimation 
of the corresponding star formation activity. Dust correc- 
tion factors are still very uncertain mainly because of the 
unknown shape of the dust UV-extinction law (Pettini et al. 
1997). The amount of star formation that is hidden by dust 
over the entire history of the universe will be constrained 



■a. 




redshift 

Figure 1. Evolution of the log of the SFR density 
(M yr- 1 Mpc -3 ) for Q = 1, (A = 0), h = 0.5 and a Salpeter 
IMF. The data points are taken from Gallego et al. (1995) (empty 
dot), Treyer et al. (1997) (filled triangle), Lilly et al. (1996) (filled 
dots), Connolly et al. (1997) (empty squares) and the HDF (filled 
squares). The solid line represents the SFR density evolution as- 
sumed in the present analysis. The fit has been kindly provided 
by P. Madau. 



by future observations of the cosmic infrared background 
(Franceschini et al. 1997; Guiderdoni et al. 1997) . So far, the 
upper limits based on the analysis of COBE-FIRAS residu- 
als poorly constrain the high-z evolution of the SFR density, 
leaving the effect of dust obscuration still as a matter of dis- 
cussion. 

The model for the global SFR evolution presented in 
Figure 1 assumes an extinction law similar to that which 
applies to stars in the Small Magellanic Cloud (SMC) and 
an amount of dust which results in an upwards correction of 
the comoving UV luminosity by a factor 1.4 at 2800A and 
2.1 at 1500A (Madau, Pozzetti & Dickinson 1997). 

An extensive investigation of this model has been ac- 
complished by Madau, Pozzetti & Dickinson (1997). For 
simplicity, the effects of cosmic chemical evolution is not 
included and all population synthesis codes assume solar 
metallicity, thus generating colors that are slightly too red 
for primeval galaxies with low metallicity. However, an en- 
couraging similar trend is being indicated by the gaseous 
and chemical evolution of the intergalactic gas as delineated 
by the studies of QSO absorption lines. The overall shape of 
the mean cosmological mass density contributed by Lyman- 
q absorbers indicates that the bulk of star formation activity 
should take place at z<2 — 3 (Lanzetta et al. 1996; Storrie- 
Lombardi et al. 1996). 

Likewise, the study of the metallic clouds in combi- 
nation with chemical evolution models leads to a volume- 
averaged star formation history in agreement with the evo- 
lution inferred from galaxy surveys (Pei & Fall 1995). 

A further important check on the inferred SFR history 
of field galaxies comes from the study of the predicted emis- 
sion properties at UV, optical and far-infrared wavelengths. 
By construction, the model produces the right amount of ul- 
traviolet light. While for a Salpeter IMF the agreement with 
the observational data still holds at longer wavelengths, a 
Scalo IMF produces too much long wavelength light by the 
present epoch. Moreover, the SFR history plotted in Figure 
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1 appears able to account for the entire background light 
recorded in the galaxy counts down to the very faint mag- 
nitudes probed by the HDF and produces visible mass-to- 
light ratios at the present epoch which are consistent with 
the values observed in nearby galaxies (Madau, Pozzetti & 
Dickinson 1997). 

Finally, a broad class of hierarchical clustering models 
predicts a SFR history in good agreement with the obser- 
vationally estimated one, giving to the model a significant 
theoretical support (Baugh et al. 1997). 



3 THE RATE OF BLACK HOLE COLLAPSES 

Using the SFR density given in Figure 1, we calculate the 
evolution of the rate of core-collapse supernovae Rsn(z), 
i.e. the number of events occurring per unit time within the 
comoving volume out to redshift z, 



Rsn(z) = f P*{z) dz' [ I 

JO dZ J M„ 



$(M) dM. 



(3) 



Here p*(z) is the SFR density, dV the comoving volume 
element and $(M) the IMF. The lower limit of the mass 
range, M p , depends on the specific nature of the supernovae 
considered. Numerical studies show that single stars with 
masses > 8Mq evolve rather smoothly through all phases 
of nuclear burning ending their life as supernovae. For this 
class, which includes Type II and Types Ibc supernovae, the 
explosion occurs through the gravitational collapse of their 
core, leaving behind a neutron star or a black hole (e.g. 
Ruiz-Lapuente 1997). The smallest progenitor mass which 
is expected to lead to a black hole ranges from 18Mq to 
30Mq, depending on the iron core masses expected from 
stellar evolution calculations (Timmes, Woosley & Weaver 
1995; Woosley & Timmes 1996). In the calculations to follow 
we will adopt a reference value of 25M© and we will discuss 
the implications of higher values of M p on our results in 
Section 6. 

As previously mentioned, p, (z) has been calculated by 
assuming a Salpeter IMF, which we will accordingly use for 
our calculations, 



$(M) oc M 



-(l+x) 



with x = 1.35, 



normalized through the relation 



M $(M) dM 



(4) 



(•») 



with Mi = O.1M and M u = 125M . 

As recently suggested by Madau (1998), the evaluation 
of the rate of core-collapse supernovae based on the SFR 
density previously discussed and on an assumed universal 
IMF is largely independent of the choice of the IMF. Using 
an extreme approach, we may assert that we do not need an 
IMF to deduce the rate of supernovae explosion. In fact, the 
SFR density is inferred from the rest-frame UV continuum 
emission, and the stars which are responsible for this emis- 
sion are the more massive ones, which, in turn, are those 
that at the end of their evolution give origin to core-collapse 
supernovae. 

The model we use for the SFR history assumes a flat 
cosmology (fin = 1) with vanishing cosmological constant 




Black Holes • 
Core-Collapse Sne 



1 2 3 4 5 

redshift 

Figure 2. Evolution of the rate of supernovae as a function of red- 
shift. The upper curve represents the rate of core-collapse super- 
novae (Type II) and assumes a lower mass cut-off of M p = 8 Mq. 
The lower curve represents only those events that lead to black 
hole collapses and assumes M p = 25 Mq . 



and h — 0.5, which will be the values of the cosmological 
parameters we will adopt in the rest of this paper. In this 
case, the comoving volume element is related to z through, 



Sr =167r (i|) 



! [1 -(!+*) 



-1/212 



(1 + z) 3 / 2 



(6) 



We have evaluated the integral (Q) by using the star for- 
mation rate density p*{z) plotted in Figure 1, and eqs. (Q) 
and (|^), for two values of M v . The lowest value of the mass 
cutoff M p = 8Mq corresponds to the total rate of core- 
collapses which produce either neutron stars or black holes. 
M p — 25Mq is a possible mass cutoff for the production of 
black holes at the present state of our knowledge (Woosley 
& Timmes 1996). The results are plotted in Figure 2. 

The predictions for the rate of supernovae are in good 
agreement with the observed local values (Cappellaro et 
al. 1997). At higher redshift, z ~ 1, the detection of Type 
II events must await future experiments, such as the Next 
Generation Space Telescope (NGST) (Madau, Delia Valle, 
Panagia 1998; Sadat et al. 1997). It should be noted that 
the rate shown in Figure 2, which is integrated over the co- 
moving volume as indicated in eq. (Eh, tends to saturate to 
a constant value: this follows from the drop of the SFR den- 
sity at high redshift. Conversely, the rate per unit comoving 
volume, would trace the rise and drop of the SFR density 
evolution. 

Finally, the total number of supernovae explosions per 
unit time leading to black hole formation is obtained by 
integrating eq. (J3J) for z = oo, 



Rbh = 4.74 events/s 



for 



M D 



25 M { 



©• 



(7) 



A significant quantity which indicates whether the collective 
effect of the bursts of gravitational waves emitted in such col- 
lapses generates a continuous background is the duty cycle. 
It is defined as 



D(z)= / dR BH (z) Ar G w (1 + 



(8) 



where D(z) is to be interpreted as the duty cycle now of all 
bursts generated back to redshift z and Atqw is the average 
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time duration of single bursts at the emission, which we will 
assume to be ~ 1 ms. The total duty cycle is 



D = 1.57 x 10" 



for 



M p = 25M Q . 



(9) 



The smallness of the duty cycle implies that our stochas- 
tic gravitational-wave background cannot be considered as 
a continuous one but rather as a shot-noise process, consist- 
ing in a sequence of widely spaced bursts with mean time 
separation of ~0.2 s much larger than the typical duration 
of each single burst, which is of order a few ms. 

It has recently been suggested that supernovae occur- 
ring soon after galaxy formation could contribute a duty cy- 
cle of order unity, creating a nearly continuous background 
(Blair & Ju 1997). These investigations are based on a super- 
nova rate obtained from its observationally estimated local 
value through an integration of the source count equation. 
Therefore, the important effect of the global SFR evolution 
was not included, leading to an overestimation of the rate. 

The implication of the evolving galaxy population has 
been investigated by Kosenko & Postnov (1998) to estimate 
the stochastic gravitational background produced by extra- 
galactic merging binary white dwarfs. They show that when 
the global SFR evolution is included, the level of this extra- 
galactic background can be comparable with the correspond- 
ing galactic background signal. 



4 THE ENERGY SPECTRUM OF A 

CORE-COLLAPSE TO A BLACK HOLE 

The energy spectra that we shall use to model the process 
of collapse to a black hole have been obtained by Stark & 
Piran (1985, 1986). They numerically integrated the fully 
non-linear system of Einstein + hydrodynamics equations 
describing the collapse of a star with a polytropic equation 
of state p = KfF , and adiabatic index T = 2. The initial 
configuration is a spherically symmetric very compact star 
with central density 1.9 x 10 15 (M core /M Q )~ 2 g cm" 3 and 
radius 8.8 x 10 J (M core /MQ) cm which, as a consequence 
of a sudden reduction to a fraction f v (with f p — 0.01 or 
0.4) of the equilibrium central pressure, starts to collapse. 
Simultaneously, the star is given an angular momentum dis- 
tribution, approximating rigid-body rotation. 

Thus, the collapse which is considered is that of a com- 
pact naked core. The studies on stellar evolution show that 
compact cores (typically degenerate iron cores) which form 
in the interior of massive stars, are surrounded by layers of 
lighter material. The dynamics of this envelope is not consid- 
ered in the simulation of Stark and Piran. In the concluding 
remarks we will return on this point, which, however, is not 
of substantial relevance for our present calculations. 

Stark and Piran find that if the total rotation parameter 
a = J I (GM 2 ore / c) exceeds a critical value 



a crlt = 1.2 ± 0.2 if f p = 0.01, 

a crlt = 0.80 ± 0.05 if f p = 0.4, 



(10) 



the rotational energy dominates, the star bounces and no 
collapse occurs. Conversely, if a < a cr i t a black hole forms, 
and a strong burst of gravitational waves is emitted, with 
the following characteristics. 

Let hgg (t,r,0, (f>) and hgj, (t,r,9, (ft) be the physical 




> 
< 



Figure 3. The average energy flux emitted during the axisym- 
metric collapse of a rotating, polytropic star to a black hole of 
M core = 1.5 Mq at a distance of 15 Mpc. The three curves cor- 
respond to assigned values of the rotational parameter and refer 
to f p = 0.01. 



components of the radiative part of the metric tensor eval- 
uated in the transverse-traceless gauge. The energy flux ra- 
diated per unit time by the collapsing source is given by the 
(Ol)-component of the stress-energy tensor of gravitational 
waves, 

dE 



Tor = 



dSdt 



-L { [h T eI{t, r, </>, 9)] 2 +[h T e J{t, r, </>, 9)] 2 } (11) 



(in geometric units, c = G = 1), where x° = t,x 1 = r, x 2 



hu, v (w, r, 



yTT /, 



' dt, 



(12) 



and apply Parseval's theorem, the energy flux per unit fre- 
quency can be defined as 



{\h^f(LU,r,(b,9)\ 2 + \hU(uj,r,(b,e)\ 2 }, (13) 



dE uj 2 
dSduj ~ 32tt 2 
and the energy flux per unit solid angle, 

sib = {v&^m? + itf?w.M)i a } • d4) 

For ui > 0, eq. ([l3|) and eq. dli| ) must be multiplied by 
a factor of 2 in order to account for the folding of negative 
frequencies into positive. We finally define an avera ge e nergy 
flux per unit frequency, f(f), by integrating eq. (|1J) over 
the solid angle, and dividing by 47rr 2 

/■2-7T 



/O) 



47rr 2 



(- 



,IB ) sm9d9. 



\ dQdu 



(15) 



'o Jo 

To express f[y) in physical units, eq. (^) must be multi- 
plied by the factor c 3 /G. 

Let us now consider the collapses studied by Stark and 
Piran that belong to the case f p = 0.01. During an axisym- 
metric collapse, the energy is emitted in both polarization 
modes but since h+ w 10 h x , we shall consider only the 
contribution of the h+ component. The function f(y) is 
plotted in Figure 3 for a black hole of M core = 1.5Mq 
which forms in the Virgo cluster, for assigned values of the 
angular momentum. The quantity f(y) attains a sharp max- 
imum for some value of the frequency which depends on the 
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angular momentum; for instance, we find i/ m „ = 8.6 kHz 
for a = 0.79 and v max = 9.2 kHz for a = 0.94. 

It is interesting to compare the location of this maxi- 
mum with the frequency of the lowest m = quasi-normal 
mode of a rotating black hole with the same mass and angu- 
lar momentum (e.g. Leaver 1986). For i — 2 we find vq = 8.6 
kHz, for a = 0.79, and vq = 9.0 kHz, for a = 0.94; for a = 
it would be uo ~ 8 kHz. Thus the peak of f(v) is located at 
a frequency which is very close to the frequency of the lowest 
m = quasi-normal mode. This means that a substantial 
fraction of energy will be emitted after the black hole has 
formed: it will oscillate in its quasi-normal modes until its 
residual mechanical energy is radiated away in gravitational 
waves. 

For values of the rotation parameter in the range 0.2 < 
a < 0.8, the amplitude of the peak of f(v) scales as a 4 . 
For high values of the angular momentum, the star becomes 
flattened into the equatorial plane and then bounces verti- 
cally, but still continues to collapse inward until the black 
hole is formed. In this case the amplitude of f(y) in the low 
frequency region increases, and a further peak appears, the 
amplitude of which may become comparable to that of the 
peak corresponding to the quasi-normal modes. 

According to Stark and Piran the main characteristics 
of the energy spectra and waveforms for f p = 0.01 and 
for f p = 0.4 are similar. If f p = 0.01 the star becomes 
more flattened and the collapse is slightly faster. The ratio 
between the maximum of the energy spectrum for f p = 0.01 
and for f p — 0.4, for angular momenta close to the critical 
value, is ~ 5. 

In general, the efficiency of the process of axisymmetric 
core-collapse to a black hole studied by Stark and Piran 
is AEGw/McoreC 2 < 7 x 10 -4 . It should be remembered 
that less symmetric situations may result in a more efficient 
production of gravitational waves. 



5 SPECTRAL PROPERTIES OF THE 
STOCHASTIC BACKGROUND 

The spectral energy density of the stochastic background 
produced by those gravitational collapses that led to black 
holes formation, can be obtained by integrating over the 
allowed range of masses and redshifts the differential rate of 
black hole formation, cLRbh{M, z), multiplied by the energy 
spectrum of a single event Thus, we need to know how the 
energy flux produced by a star of mass M and angular 
momentum a, which collapses at a given redshift z, would 
be observed today. 

The average energy flux per unit frequency in the source 
rest frame, which is plotted in Figure 3 for a core of 1.5 Mq 
at the distance of 15 Mpc can be written as follows 



f{u) = A{a,v) 



Mr, 



1.5M e 



d 



15Mpc 



(16) 



where d is the distance from the source in the source rest 
frame, and 



t The smallness of the total duty cycle implies that, on average, 
different bursts do not overlap, so that interference among the 
waves generated by different collapses is unlikely to occur. 



'emission 



JMcor 



1.5M 

Such a signal, if observed here and now, would give 



f(v bs) = A(a,v) 



( d L {z) V 
1.5M Qy / ^15Mpcy 



(17) 



(18) 



where the observational frequency v t, s is related to v 
through 



v 1.5M, 



Vobs = 



and 







(l + z) M cor 



2c 



(19) 



(20) 



4(^j(i+ 2 )[i-(i+ z r 1/2 i 
Ha 

is the luminosity distance. In the following we shall assume 
that the fraction of the progenitor star which collapses is 
M core = a M with a = 0.1. 

The differential rate of black hole collapse is given by 
[see eq. (§)] 



dR BH {M,z) = p*(z) $(M) dMdz. 
Thus the spectral energy density is 

coo fM u 

f(v obs )dR B H(M,z) 



dE 



dtdSdv 



(21) 



(22) 



Since we do not know the distribution of angular momenta of 
the formed black holes, we evaluate the spectral energy den- 
sity for assigned values of the rotation parameter. We con- 
sider three possible values of a, namely a — 0.42, 0.79, 0.94, 
and give the corresponding functions A(a, v) as an input to 
the numerical code. For each chosen value of a, we solve the 
integral (^) as a function of the observation frequency v a i a 

The statistical description of a stochastic gravitational 
background relies on three basic assumptions, i.e. that the 
signal is isotropic, stationary and unpolarized (Allen & 
Romano 1998, Maggiore 1998). In our present case, the 
background is produced by a large number of indepen- 
dent gravity-wave sources extending over the redshift range 
< z < 4 — 5. Thus, the first two assumptions are cer- 
tainly satisfied and, assuming that the sources have random 
orientations, the radiation incident on the detector has sta- 
tistically equivalent polarization components (Thorne 1987). 
From the spectral energy density we shall evaluate the cor- 
responding values of the spectral strain amplitude, 



VSn(a^o b s)= (f|4-) ( 



dE \ 



1/2 



, dtdSdv 



) 



(23) 



and of the closure energy density of gravitational waves, 



flGw{a, Vobs 



Voba dE 4 7T 2 



u 3 obs S h (a,u). (24) 



C Per 

dtdSdv 3 Hi 
The results are plotted in Figures 4, 5 and 6 



6 DEPENDENCE ON THE ASSUMED MODEL 
PARAMETERS 

Let us now discuss how the spectral energy density and the 
spectral strain amplitude computed in Section 5 depend on 
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Figure 4. The spectral energy density (dE / dSdtdu) is plotted 
as a function of the observational frequency. The three curves 
correspond to assigned values of the rotational parameter. 




1/2 

Figure 5. The spectral strain amplitude SJ corresponding 
to the spectral energy density given in Figure 4 is plotted as a 
function of the observational frequency. 




Figure 6. The function £Igw corresponding to the spectral en- 
ergy density plotted in Figure 4. 



some of the parameters we use in our calculations. First of 
all the mass cutoffs. 

The lower cutoff on the mass of the progenitor star 
whose collapse produces a black hole, and which we set at 
M p = 25Mq, depends on the stellar evolution model upon 
which the calculations describing the life of massive stars are 
based, as well as on the effect of fallback during the explo- 
sion (Woosley & Timmes 1996 and reference therein). Thus, 
we have investigated the consequences that a higher value of 
this cutoff, namely M v = 30 Mq , would introduce in our re- 
sults and we have similarly considered the effect of reducing 
the dynamical range lowering the upper mass cutoff from 
125 Mq to 60 Mq . The various dynamical ranges considered 
are listed in Table 1, where we also give the values of the rate 
of collapses to a black hole and the duty cycle. We see that 
if the lower mass cutoff is raised to M v = 30 Mq the values 
of the total rate and duty cycle are a factor ~ 1.33 lower 
than in the M v = 25 Mq case. The effect of lowering the up- 
per mass cutoff to 60 Mq is less significant. Infact, although 
the percentage of progenitors having mass between 60 and 
125A40 is about 38%, when we compute the rate as in eq. 
([]), since the IMF goes like ~ M -2 ' 35 , bigger masses will 
give smaller contributions to the integral, and the change in 
rate and duty cycle will be small. 

As for the spectral energy density, since the frequen- 
cies of the energy spectrum emitted during a single collapse 
are inversely proportional to the mass of the collapsing core 
(which we have assumed to be roughly 10 % of the progen- 
itor mass), the major effect of shifting the dynamical range 
towards higher values of the progenitor mass will be to re- 
duce the power at high frequency and decrease the overall 
amplitude. Similarly, the effect of reducing the upper mass 
cutoff is that of loosing some power at low frequencies. In 
addition, from Figure 7 we see that the effect of increasing 
M p to 30 Mq is less significant than that of lowering M u to 
60 Mq. The reason is immediately understood if we look at 
eq. (G2f) from which dE / dtdSdv is computed. It involves an 
integral over the mass which, due to the explicit dependence 
of f(v obs ) and of $(M) on M, goes like M" 35 . 

As for the spectral strain amplitude, the l/u b s factor 
in eq. (^3|) introduces a factor M, through eq. (|l^), which 
makes negligible the change in ySh when the lower mass 
cutoff is raised to 30 Mq . Conversely, the effect of changing 
the upper mass cutoff is to decrease the amplitude of the 
peak by almost a factor of 2. 

The results of this discussion are plotted in Figures 7 
and 8 for a single value of the rotational parameter a — 

0. 79 and the for various combinations of the mass ranges 
given in Table 1. In Table 2 we give the frequency at which 
the maximum of [dE / dtdSdv) and \/~Sh is located in the 
various cases, and the corresponding amplitudes. 

It is also interesting to see how the strain amplitude 
changes if we choose different values of the parameter a, 

1. e. of the fraction of the progenitor star which collapses. 
If a increases, the mass of the collapsing core increases. As 
a consequence, the emitted energy, which is proportional of 
M 2 ore , increases, and the frequency of the maximum, which 
is inversely proportional to M core decreases. This results in 

1/2 

an overall higher amplitude of S h , and in the moving of 
the peak frequency toward smaller values. The results are 
shown in Figure 9. 

We have finally analysed the consequences of a differ- 
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Table 1. The different combinations of low and high mass cutoffs considered in the analysis of the dependence of the results on the 
assumed model parameters. 



Model 


M P /M Q 


M„/M 


Rbh (events/s) 


D 


case A 


25 


125 


4.74 


1.57 x 1(T 2 


case B 


30 


125 


3.57 


1.19 x 1(T 2 


case C 


25 


60 


3.71 


1.23 x icr 2 



1/2 

Table 2. The values of the frequencies at which dE/dtdSdu and S h have a maximum, and the corresponding amplitude (respectively 
in erg cm - 2 s- 1 Hz" 1 and in Hz -1 / 2 ) are tabulated for different combinations of the upper and lower cutoff parameters. 

Model M P /M Q M up /M Q iw(Hz) (3^;)^ w(Hz) v<S h max (Hz" 1 / 2 ) 



case A 25 125 1.05 x 10 3 

case B 30 125 1.05 x 10 3 

case C 25 60 2.55 x 10 3 



5.0 x lO" 11 3.12 x 10 2 5.3 x 10" 28 
4.9 x lO" 11 3.12 x 10 2 5.4 x 10" 28 
4.0 x KT 11 7.56 x 10 2 2.4 x 10~ 28 





Figure 7. The spectral energy density dE / dSdtdu is plotted 
for various combinations of the mass cutoffs (see Table 1). The 
curves correspond to a unique value of the rotational parameter 
a = 0.79. 



1/2 

Figure 9. The spectral strain amplitude S h is plotted for three 
different values of the parameter a which defines the fraction of 
the progenitor star which collapses. The curves correspond to 
a = 0.79. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
v (kHz) 

1/2 

Figure 8. The spectral strain amplitude SJ corresponding to 
the spectral energy density given in Figure 6. 



ent evolution history of the SFR density at high redsliift, 
assuming a constant tail at z > 1 — 2 at approximately the 
amplitude of the maximum. This evolutionary behaviour is 
still consistent with the upper limits found in the analysis of 
COBE-FIRAS residuals (once the black body contribution 
of the cosmic microwave background has been subtracted 
off) in the infrared and is predicted by models of galaxy 
evolution in which dust obscuration significantly attenuates 
the UV emission at high redshifts (Franceschini et al. 1997). 
The overall effect of an increased SFR at high redshifts is 
negligible on the spectral energy density dE / dSdtdv as well 
as on Q.gw and sligthly increases the amplitude of ■s/Sh at 
low frequencies (see Figure 10). 

Thus, the results of our analysis are not seriously af- 
fected by the uncertain amount of dust extinction, and the 
possibility that UV-optical data may lead to an underesti- 
mation of the global SFR at high redshifts does not inval- 
idate our main conclusions. The reason for the negligible 
role played by high redshift sources is that the average en- 
ergy flux of gravitational waves contributed by each source 
decreases as the inverse of the squared luminosity distance. 
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X 4 



UV-optical SFR density 
constant SFR density at z>l-2 





1 /2 

Figure 10. The spectral strain amplitude S h is plotted for two 
different evolutionary models for the SFR density: the continuous 
line refers to the SFR shown in Figure 1, the dotted line uses the 
same SFR but with a constant tail at z > 2. (see text). The 
curves correspond to a unique value of the rotational parameter 
a = 0.79. 



7 STATISTICAL PROPERTIES OF THE 
STOCHASTIC BACKGROUND 

In Sections 5 and 6, the stochastic background produced by 
gravitational collapses to black holes throughout the uni- 
verse was considered to be continuous. This is indeed the 
case as long as an observational run is much longer than the 
mean temporal interval between two successive bursts. This 
is given by the inverse of the total rate of events and can 
be estimated to be ~ 0.2 s (see Section 3). Likewise primor- 
dial stochastic backgrounds, the statistical properties of the 
signal will be similar to the detector noise but with lower 
amplitude and the optimal detection strategy will consist 
in performing a correlation between at least two detectors 
(Allen 1996; Allen & Romano 1997). 

As an alternative strategy, one could try to exploit the 
characteristic shot-noise structure of the signal and design a 
specific algorithm to extract it from a single detector noise. 
We plan to discuss this possibility in a forthcoming paper. 
As a preliminary step, we shall now discuss some statistical 
properties of this background. We have computed the num- 
ber of individual signals that are expected to have a max- 
imum amplitude hmax above a given threshold, h t , within 
an observational run AT a t s 



J^AT obs (h m ax > h t ) = V(hmax > ht) AT obs R B H, 



(25) 



where V(h ma x > ht) is the probability that h max is above 
the threshold, and Rbh is the total rate of collapses to a 
black hole. 

Following Stark and Piran, the expression of the maxi- 
mum amplitude of the gravitational wave emitted in a single 
black hole collapse can be written as 



hma X {M,z) = C(a) 



aM 



( d L (z) \ 1 



1.5 M Q yi5Mpc J 



(26) 



where C(a) is a constant term which is related to the rota- 
tional parameter a through, 



C{a) 



4.79 x 10" 
2.86 x 10" 



if a < 0.8 
otherwise. 



Figure 11. The probability of detecting a signal with maximum 
observable strain amplitude h max greater than a threshold value 
ht, is plotted as a function of ht for assigned values of the rota- 
tional parameter. 



Thus, h max (M,z) depends on the mass of the progenitor 
star and on the redshift at which the collapse occurs. There- 
fore, its probability distribution can be obtained from the 
probability distributions of the mass and redshift p(M, z) 
by performing the following integral, 



p(h) 



dMdz p(M,z) 8 D (h max (M,z) - h), 



(27) 



where h is any value of the random variable h max and 8d is 
the Dirac delta function. 

The function p(M, z) is computed in the following way. 
The number of black hole collapses with progenitor mass in 
the range (M, M + dM) occurring at a redshift between 
(z, z + dz) is, 



dM = $(M) dM V(z) — dz. 

dz 



(28) 



Here ^(z) is the mass of gas that goes into stars per unit 
comoving volume, i.e., 



*(*) = / P.(*) 



dt 



dz' 



dz , 



(29) 



where, in our case the cosmic time is related to z through 
dt/dz = -Hq 1 (1 + z)~ 5/2 . Thus, the probability that a 
black hole collapse occurs with progenitor mass in the range 
(M, M + dM) and redshift in the range (z, z + dz) is, 

p(M, z) dM dz = (30) 

J\tot 

where Mtot is the total number of events. The probability 
of a signal with h max > h t , 



V(hmax > h t ) 



dh p(h) 



(31) 



is plotted in Figure 11 as a function of ht for the same values 
of the rotational parameter considered in Section 5. 

In Figures 12 and 13 we show the corresponding number 
of events expected A/at d6s (hmax > ht) for an observational 
run of 1 s and 1 yr respectively. 
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Figure 12. The number of events with hmax > ht in an obser- 
vational run of 1 sec as a function of ht for three assigned values 
of the rotational parameter. 




Figure 13. The number of events with hmax > ht in an obser- 
vational run of 1 yr as a function of ht for three assigned values 
of the rotational parameter. 

8 CONCLUDING REMARKS 

From the analysis carried out in this paper it emerges that 
a reliable estimate of the stochastic background of gravita- 
tional waves produced by a cosmological population of as- 
trophysical sources cannot set the important effect of SFR 
evolution aside. 

The extraordinary advances attained in observational 
cosmology are leading to a coherent picture for the global 
star formation history in field galaxies that can be used to 
infer the redshift evolution of the rate of gravitational wave 
sources. 

More importantly, the resulting gravitational signal ap- 
pears to be insensitive to the uncertainties that still trouble 
the high redshift SFR observations (as a consequence of the 
uncertain amount of dust obscuration), as the main contri- 
bution to the resulting Qgw and V~Sh comes from low-to- 
intermediate redshift sources. 

We have restricted our investigation to a specific cos- 
mological background model with Qo = 1, A = 0, Ho = 
50 kms~ 1 Mpc -3 . We plan to extend the analysis to differ- 
ent choices of background cosmology in a subsequent paper. 
We can anticipate that, if we assume that the cosmological 



constant is zero, the SFR density decreases with Qo- At the 
same time, the geometrical effect on the comoving volume el- 
ement increases the rate of black hole collapses compared to 
a flat background, although the decrease of the single event 
spectra with the square of the inverse luminosity distance 
partly cancels the effect. Altogether we would then expect a 
modest increase of amplitude of the resulting gravity-wave 
background. 

The spectral properties of the background of gravita- 
tional waves discussed in this paper depend on the energy 
spectrum we have chosen as representative of the process 
of gravitational collapse to a black hole. As mentioned in 
Section 4, it was derived by Stark & Piran (1985, 1986) and 
it refers to a simplified model of star, namely a naked com- 
pact core with a polytropic equation of state. At present, 
detailed information on the gravitational radiation emitted 
in collapses of stars with more realistic equations of state, 
and derived by a fully relativistic approach as in the Stark- 
Piran model, are not available. 

As far as the "naked core" assumption is concerned, it 
should be mentioned that detailed studies on the structure 
of stars that evolve into Type II supernovae, show that quite 
a general pre-collapse configuration is composed by a degen- 
erate iron core surrounded by layers of matter, remnants of 
the earlier nuclear burning processes (Woosley & Weaver 
1995) However, these layers are not expected to contribute 
significantly to the emitted radiation, because they do not 
participate in the core bounce (see for example the discus- 
sion in Section VII of Seidel 1991). 

The preminent features of the Stark-Piran energy spec- 
trum are (see Figure 3): 

• The presence of a peak at the frequency of the quasi- 
normal modes of the formed black hole. 

• The dependence of the efficiency on the fourth power 
of the angular momentum of the collapsing star. 

• The presence of a further pronounced peak at a lower 
frequency, if the angular momentum is sufficiently high. 

It is interesting to ask to what extent are these features 
general, and how do they reflect in the structure of the 
quantities we compute, for instance the strain amplitude, 
which is of direct interest for gravitational detectors. Grav- 
itational signals emitted during the collapse of naked cores 
to a black hole have been studied in several papers with a 
perturbative approach: a spherically symmetric collapsing 
star is perturbed axisymmetrically, and the emitted radi- 
ation is computed by solving the linearized equations (see 
Ferrari & Palomba 1998, for a recent review). The amount 
of gravitational energy estimated by the perturbative ap- 
proach is obviously lower than that predicted by Stark and 
Piran. However, the presence of the peak of the quasi-normal 
modes can be considered as a general feature, unless the 
dynamics of the fluid is dominated by very strong internal 
pressure or by very long bounces, which basically slow down 
the collapse to such an extent that the formed black hole has 
no residual mechanical energy to radiate away in its quasi- 
normal modes (Seidel 1991). The "two-peak" structure of 
the energy spectrum can be mantained also if the star is 
non-rotating (Seidel 1991), if bounces occur during the col- 
lapse. The dependence of the emitted energy on the fourth 
power of the angular momentum is found also in the per- 
turbative approach (Cunningham, Price & Moncrief 1980). 
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In all cases the frequency scales as the inverse of the black 
hole mass, and the energy radiated as the second power of 
the mass. Thus, the features of the energy spectrum we use 
to model each single event are likely to reasonably represent 
a generic situation. 

The strain amplitude we computed from the Stark- 
Piran energy spectrum is shown in Figure 5. It exhibits a 
peak at some frequency which is the reminiscent of the peak 
of the quasi-normal modes of the formed black holes. Since 
we do not know the angular momentum distribution of the 
black holes, each curve in that figure is obtained by assum- 
ing that all black holes that form have the same value of the 
rotation parameter, and the location of the peak of \fSh 
depends on that value. We can summarize the above discus- 
sion by saying that the strain amplitude associated to the 
background under consideration should present a peak at a 
frequency between ~ 230 and 340 Hz, whereas the expected 
maximum amplitude should be of order 1 — 8 x 10 -28 . 
However, since our results are based on the assumption that 
each collapse takes place axisymmetrically the maximum 
strain amplitude may be underestimated: non-axisymmetric 
collapse would be more efficient in producing gravitational 
waves. The strain amplitude plotted in Figure 5 was derived 
for the following set of parameters: a — 0.1, M p = 25Mq 
and M u — 125M©. The dependence on different choices for 
the progenitor mass range (M p , M u ) was considered in sec- 
tion 6. The results are shown in Figure 8 and in Table 2. 
For a lower value of the upper mass cutoff, the frequency 
of the maximum is shifted by a factor ~ 2.4 towards higher 
frequencies whereas the maximum amplitude is lowered by 
a factor ~ 2.2. Similarly, Figure 9 shows the effects of vari- 
ating the fraction of the progenitor mass which participates 
to the collapse (a). We can consider these curves as an in- 
dication of the uncertainties introduced by the lack of con- 
sistent results in the literature. However, reasonable values 
of the parameters a, M p and M u should not be too far from 
the ranges we have adopted. Thus, we may conclude that, 
according to our study, the spectral strain amplitude pro- 
duced by the background of supernova explosions leading 
to black hole formation should have a maximum at a few 
hundred Hz with an amplitude ranging between 10 -28 and 
10~ 27 Hz _1//2 . The corresponding closure density, Qgw, has 
a maximum amplitude ranging between 10" 11 and 10" 10 
in the interval ~ 1.5 — 2.5 kHz. This is beyond the sensi- 
tivity that can be obtained by cross-correlating the outputs 
of two gravity-wave detectors of the first generation experi- 
ments (VIRGO- VIRGO, LIGO-LIGO). However, for a con- 
stant frequency spectrum the advanced LIGO detector pair 
aims at Q,gw ~2x 10" 10 (Allen & Romano 1998), a level 
comparable to our signal. 

Finally, we have made a preliminary statistical analysis 
of our background, e.g. in terms of the probability that the 
signal is above a given threshold. More details will be consid- 
ered in a forthcoming paper. Let us stress, once again, that 
the smallness of the duty cycle, which makes our stochas- 
tic background a non-continuous one, implies a non-trivial 
statistics of the signal arrival times which will be a distinc- 
tive feature of the process that can be in principle exploited 
to design optimal detection strategies. 
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